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O Abstract: 

(— I The spectral statistics in the strongly chaotic cardioid billiard are studied. The analysis is 
O based on the first 11000 quantal energy levels for odd and even symmetry respectively. It is 
found that the level-spacing distribution is in good agreement with the GOE distribution of 
random-matrix theory. In case of the number variance and rigidity we observe agreement with 
the random-matrix model for short-range correlations only, whereas for long-range correlations 
both statistics saturate in agreement with semiclassical expectations. Furthermore the con- 
jecture that for classically chaotic systems the normalized mode fiuctuations have a universal 
Gaussian distribution with unit variance is tested and found to be in very good agreement 
for both symmetry classes. By means of the Gutzwiller trace formula the trace of the cosine- 
modulated heat kernel is studied. Since the billiard boundary is focusing there are conjugate 
points giving rise to zeros at the locations of the periodic orbits instead of exclusively Gaussian 
peaks. 
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1 Introduction 



Today the important role played by chaos in nonlinear dynamical systems is generally appreci- 
ated. The most striking property of deterministic chaos is the sensitive dependence on initial 
conditions such that neighbouring trajectories in phase space separate at an exponential rate. 
As a result the long-time behaviour of a strongly chaotic system is unpredictable. There the 
fundamental question arises whether this well-established phenomenon of classical chaos has 
an analogue in the quantum world that could be called "quantum chaos". By this the follow- 
ing is meant: Given a dynamical system which is strongly chaotic, i.e., ergodic, mixing and 
a K-system, is there any manifestation in the corresponding quantal system which betrays its 
chaotic character? (For an authoritative review, see ref. [1].) If one were to identify unique 
fingerprints of classical chaos in quantum mechanics, one could use these to define quantum 
chaos. Ideally, classically chaotic systems should be characterized by a random behaviour of 
these fingerprints that qualify the system to be called "chaotic" also in quantum mechanics. 

For bound conservative systems, the quantum mechanical time evolution is almost periodic, 
in the sense of Harald Bohr's theory of almost periodic functions, due to the discrete spectrum 
of the time-evolution operator. One thus observes no sensitive dependence on initial conditions 
in the long-time behaviour of quantum mechanics. This is in contrast to classical systems whose 
time evolution is ruled by the Liouville operator. If the classical dynamics is strongly chaotic, 
the spectrum of the Liouvillian has a continuous spectrum on the unit circle which leads to 
a decay of time correlations of classical observables refiecting a complete loss of information 
on the system. This fundamental difference between classical and quantum mechanics has led 
to the common belief that there hardly exists any phenomenon in quantum mechanics which 
justifies the notion of quantum chaos. 

However, instead of concentrating on the long-time behaviour, one can consider the extreme 
limit t = oo in quantum mechanics and thus study properties of stationary states, that is of 
eigenvalues and eigenfunctions of the corresponding time-independent Hamiltonian. The idea 
is that the statistical properties of the energy-level fiuctuations ("spectral statistics") and wave 
functions of a given quantum system are already determined by its classical limit, depending 
only upon whether this limit is chaotic or not. (In this paper we shall consider the spectral 
statistics only; for a discussion of the statistical properties of wave functions in chaotic systems, 
see refs. [1, 2, 3, 4, 5, 6].) It has been conjectured that the statistical properties of quantum 
energy spectra of classically integrable systems can be described [7] by Poissonian random 
processes, whereas the spectral statistics of strongly chaotic systems can be described [8] by 
the universal laws of random-matrix theory (RMT), originally proposed by Wigner, and Landau 
and Smorodinsky, and fully developed by Dyson for a better understanding of the resonances 
of compound nuclei. (See ref. [9] for a collection of the original papers, and refs. [1, 10, 11] 
for recent reviews.) The random-matrix model for spectral statistics has to be viewed as a 
purely phenomenological one, in contrast to random-matrix theory which constitutes an exact 
mathematical theory. So far there does not exist a complete theory for the spectral statistics of 
chaotic systems, nor even for integrable ones. Based on Gutzwiller's trace formula [12, 1], Berry 
carried out a semiclassical analysis [13, 14, 15] of the level fiuctuations and obtained a saturation 
of the two-point statistics for long-range correlations in contradiction with the predictions based 
on RMT; non-universal long-range correlations occur due to the non-universal behaviour of 
short periodic orbits. A recent analysis [16, 17] based on the exact Selberg trace formula has 
clearly confirmed the non-universal saturation up to correlation lengths as large as L = 700 (in 
units of the average level spacing). In addition, there exists a class of strongly chaotic systems 
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for which the spectral fluctuations nearly behave as it is expected for classically integrable 
systems. This phenomenon occurs for geodesic flows on hyperbolic manifolds (for instance, 
Riemannian surfaces with constant negative Gaussian curvature), whose fundamental groups 
are of an arithmetical origin; thus the notion of arithmetical chaos was introduced [18, 19]. The 
dynamical systems possessing arithmetical chaos violate universality in energy-level statistics 
even in the short-range regime. It thus appears that the random-matrix model does not provide 
a universal signature of classical chaos in quantum mechanics. 

Recently, a novel quantity to measure quantum chaos has been proposed and a conjecture 
about its statistical behaviour has been put forward [5, 6]. According to this conjecture there are 
unique fluctuation properties in quantum mechanics which are universal and, in a well-deflned 
sense, maximally random if the corresponding classical system is strongly chaotic. Numerical 
as well as theoretical evidence has been provided in favour of the conjecture [6, 17]. A rigorous 
proof of the conjecture would give us a clear-cut deflnition of quantum chaos in spectra. 

In order to shed more light on the relationship between classically chaotic systems and the 
corresponding quantum systems it is important to have a large number of systems for which 
it is possible to carry out extensive numerical computations of the relevant spectral statistics. 
The simplest nonlinear dynamical systems one can study are billiard systems. They consist of 
a point particle moving freely inside a curved boundary with elastic reflections at the boundary. 
These systems are specially suited for the study of a possible manifestation of classical chaos 
in quantum mechanics, because a lot of mathematical results are available. Furthermore the 
integration of the equations of motion is trivial, i.e., the geodesies are straight lines, which by no 
means does imply that it is trivial to determine all periodic orbits having periods below a given 
value. Moreover billiard systems belong to the so-called scaling systems with the semiclassical 
limit of Planck's constant /i — )• being equivalent to the high energy limit E ^ oo. 

In this paper we are concerned with a billiard system bounded by the cardioid [20]. This 
system is of special interest because it provides one of the few examples for which it is rigorously 
proven that it is strongly chaotic [21, 22, 23, 24]. Famous other examples are the Sinai billiard 
[25, 26], which is very important for the foundations of statistical mechanics, and Bunimovich's 
stadium billiard [27, 28, 29, 30]. However, the last two systems possess some non-generic 
features, because they have a family of stable periodic orbits, the so-called bouncing-ball modes. 
The cardioid billiard has no such family, and no whispering gallery orbits exist. Furthermore 
the boundary is focusing and the geometric properties are such that caustics exist. Their 
important influence will become clear in sec. 4 in the discussion of the cosine-modulated heat 
kernel. 

The paper is organized as follows. After deflning the cardioid billiard and its quantum 
mechanical version in sec. 2, we present in sections 3.1 to 3.5 a detailed analysis of the spectral 
statistics, i.e., the spectral staircase function, the (5„-statistic, the level-spacing distribution 
P(s), the number variance TP'{L) and the rigidity ^^{L). In sec. 3.6 the distribution of the 
mode-fluctuation number is calculated and found to be in excellent agreement with the recent 
conjecture [5, 6] on quantum chaos in spectra. In order to keep the number of flgures limited, 
the various spectral measures are displayed in each case for one symmetry class only, but we 
show the plots for even and odd symmetry in an alternating order. In section 4 we consider 
the trace of the cosine-modulated heat kernel additional test of the accuracy of the 

eigenvalues and as a flrst application of Gutzwiller's periodic-orbit theory. Section 5 contains 
a summary and discussion of our results. 
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2 The Cardioid Billiard 



The billiard system to be studied is given by the free motion of a point particle inside a two- 
dimensional Euclidean domain bounded by the cardioid (see figure 1) with elastic refiections 
at the boundary dft. The cardioid billiard is the limit case of a family of billiards which was 
introduced by Robnik [20]. The boundary dft of the billiards is defined by a quadratic conformal 
mapping of the unit circle 

w = u + iv = z + \z^ ■ z = e''^ , (^G[0,27r[ , Ag[0,|[ • (1) 

Starting with A = one gets a continuous deformation of the circle. The cardioid billiard is 
obtained for A = |; for this value the mapping is no longer conformal, since ^ = 1 + 2; = 
for 2; = — 1 ((^ = tt), where the cardioid has a cusp. It was proven that the cardioid billiard 
is ergodic, mixing, and a K-system [21, 22, 23]. In fact it is even a Bernoulli system, which 
follows from a recent result [24]. 




Figure 1: Full and desymmetrized cardioid billiard. D denotes Dirichlet and N Neumann 
boundary conditions. 

The corresponding quantum mechanical system is governed by the Schrodinger equation (in 
natural units fi = 2m = 1) 

-A*„(^) = En^niq) , 9 G , 

with Dirichlet condition ^„((f)=0 at the boundary dft. A is the two-dimensional Laplacian. 
Furthermore, we have the orthonormality relation 

Notice that the eigenfunctions ^n{q) can be chosen as real. For being compact the energy 
spectrum {En} is purely discrete. Since the system is invariant under refiection at the symmetry 
line [v = 0) one can classify the wave functions ^n{q) as odd and even eigenfunctions satisfying 
Dirichlet or Neumann boundary conditions respectively at the symmetry line. Therefore we 
will consider the two desymmetrized versions of the cardioid billiard only, see figure 1. 

For the following discussion it is important to keep in mind that the semiclassical limit 
/i — )• corresponds to the high energy limit En — ^ 00. 
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3 Spectral Statistics 

3.1 Spectral Staircase and Weyl's Law 

One can easily check the asymptotic behaviour of the spectrum by calculating the spectral 
staircase function (integrated level density) 

N{E) = #{n\E^ < E} (2) 

which counts the number of energy levels below E. N[E) can be divided into a smooth and an 
oscillatory part 

N{E) = N{E) + N,,,{E) . 

The mean behaviour of N[E) is asymptotically for E ^ oo described by the generalized Weyl's 
law including perimeter, corner and curvature corrections [31]. For odd ( — ) and even (+) 
symmetry respectively, one obtains 

N-(E) = —E-—VE + — , (3) 
^ ^ 16 27r 16 ' ^ ' 

N^E) = ^E-^Ve , (4) 

16 ZTT 

where the inward pointing cusp contributes like an edge with angle tt in the desymmetrized 
billiard. 




Figure 2: Spectral staircase N (E) for odd symmetry and generalized Weyl's law, equation (3). 
The insets show a magnification of the energy intervals [0,100] and [59150,59270] respectively. 

In figure 2 the spectral staircase N~[E) is shown for the first 11000 energy levels for odd 
symmetry and compared with the generalized Weyl's law N~[E)^ equation (3). Due to the large 
energy interval no difference between both curves is visible in this figure on the whole range. 
Therefore two insets at the lowest and highest end are shown. The eigenvalues were computed 
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by Prosen and Robnik [32] by using the conformal mapping diagonahzation technique which 
was introduced by Robnik [33, 34]. In ref. [34] the energy-level statistics were investigated in 
detail for the family of billiards (1) for various parameters A < | for which a rigorous proof of 
strong chaos is lacking. 

Despite the fact that N~[E) is an asymptotic law, it gives the correct mean behaviour 
down to the ground state, as it was observed in many other systems before. Figure 2 strongly 
indicates that the computed energy spectrum for odd symmetry is complete for the first 11000 
levels. In the same way we have also checked the energy spectrum for even symmetry being 
complete for the first 11000 levels. Furthermore no degeneracies of eigenvalues in the spectra 
for even and odd symmetry were found. 

In order to compare the eigenvalues of the cardioid billiard with those of other systems, it 
is necessary to unfold the spectrum by means of the generalized Weyl's law (3) or (4) 

In the following we shall analyze the unfolded spectra {E'^} for the two desymmetrized billiards, 
but the prime and superscript it will be omitted. After the process of unfolding the spectrum 
has a mean level spacing of unity, and different systems differ in the oscillating part Nosc{E) 
only. 

3.2 6„-Statistic 

A more refined way of testing whether the spectrum is complete can be performed by considering 
the fiuctuating part Nosc{En) of the spectral staircase function evaluated^ at the unfolded 
eigenvalues En 

6n := N,sc{En) := N{En) - N{En) = n - \ - E^ . (5) 

As one can see in figure 3, where the even case is shown, 8n is oscillating around zero, 
showing that the levels indeed obey the mean behaviour as described by Weyl's law. If a level 
was missing, 8n would oscillate around —1 from then on. This indicates that the spectrum for 
even symmetry is complete for the first 11000 levels. The same holds for the odd case. 

Remarkably, the (5„-statistic is even sensitive to the constant term in Weyl's law, equations 
(3) and (4), as one can infer from the mean values of 8n for which we obtain ((5) = —0.00072 for 
odd, and ((5) = —0.00217 for even symmetry. This is in the odd case more than a factor 200 less 
than the constant Codd — — 0.1875 in equation (3). In case of the constant being neglected 
or unknown, the mean value of 8n would provide a good approximation to the constant term 
in Weyl's law. 

Although the fiuctuations 8n have mean zero, a look at figure 3 suggests that their variance 
slowly increases with increasing energy. A quantitative discussion of this important observation 
will be postponed until section 3.5. 

^Here the value of the spectral staircase function (2) at an unfolded eigenvalue £"„ is defined as 



6 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 11000 

n 



Figure 3: Plot of the fluctuating part 6,, 



-^osc ( -^n ) • 



= n — 



I — En for even symmetry for the 



flrst 11000 eigenvalues 

3.3 Level-Spacing Distribution 

An important statistic measuring short-range correlations of the spectrum is the nearest- 
neighbour level-spacing statistic P{s). P[s)ds is the probability of flnding an arbitrary pair 
of nearest neighbours of energy levels with spacing s„ = En+i — En in the interval -|- ds]. 
Thus P{s) measures fluctuations of the distances between two nearest-neighbour levels. 

The level-spacing distribution for typical integrable systems is found to obey the Poissonian 
distribution P^°^^(s) = e~'% which is in clear contrast to the level repulsion observed in generic 
chaotic systems. The random-matrix model [10] predicts that for systems with time-reversal 
symmetry the level-spacing distribution P{s) is given by the distribution of the Gaussian or- 
thogonal ensemble (GOE) which is well described by a Wigner distribution 



whereas systems without such symmetry should obey the distribution of the Gaussian unitary 
ensemble (GUE) which is in good approximation given by 



Since the cardioid billiard possesses time-reversal symmetry the level-spacing distribution is 
expected to be described by a Wigner distribution. 

In flgure 4 the level-spacing distribution for odd symmetry is shown. The agreement with 
pGOE^^-j is very good. A similar result holds in the even symmetry case. 

The level-spacing distribution has the drawback that there is a loss of information due to 
the binning. Therefore it is more signiflcant to look at the cumulative level-spacing distribution 





Numerically this quantity is calculated by sorting the set of nearest-neighbour spacings 
and determining the fraction of level spacings below a given s. 





s 
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Figure 4: Histogram of the nearest-neighbour level-spacing distribution P{s) for odd symmetry. 
The dashed line is the expected Wigner distribution, the dashed-dotted line is the GUE distri- 
bution, and the dotted line shows the Poissonian distribution expected for integrable systems. 




12 3 4 

S 

Figure 5: Cumulative level spacing for even symmetry. Full line: calculated from the energy 
levels, dashed line: GOF expectation /'^^^(s), dashed-dotted line: GUF expectation /'^^^(s), 
and dotted line the Poissonian distribution 

As one can see in figure 5, where the even case is shown, the agreement between the 
cumulative level-spacing distribution and the cumulative GOF expectation is very good. The 
Kolmogorov-Smirnov test gives a significance level of V = 24% for the odd, and V = 67% for 
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the even case assuming a GOE distribution, whereas for the GUE distribution an extremely 
small significance level of < 10~^^ is obtained for both symmetry classes. 

3.4 Number Variance 

In contrast to the level-spacing statistic, the number variance allows one to investigate also 
medium- and long-range correlations of the spectrum. It is defined by 

2\ 



^\L,E):={{n{L,E)-Ly) , L >0 



E,5 

which is the local variance of the number n{L, E) = N{E + f ) - N{E - f ) of unfolded energy 
levels En in the interval [E — ^ ^ E -\- ^]. The brackets (. . .)^ ^ denote a local average with center 

E and effective width 6. Different averaging procedures are possible; here we have chosen a 
rectangular averaging 

E+S/2 

{f{E))E,s = \ j m dE . (6) 

E-512 



In random-matrix theory one obtains in the GOE case 



S^eW = ^ log(27rZ) + 7 + I + -Si2(7rZ)--Si(7rZ) 



cos(27riv) - Ci(27riv) + tt^L 



-Si{27rL] 

TT 



and for GUE 



Sgue(^) = ^ \ log(27rZ) + 7 + I - cos(27rZ) - Ci(27rZ) + tt'^L 



1 - -Si{27rL] 

TT 



(7) 



(8) 



where Si{x) and Ci(a;) are the sine and cosine integral respectively, and 7 = 0.5772 ... is Euler's 
constant. Eor a random Poissonian process one has 

^pois(-^) = L , 

which is in agreement with the general small-Z behaviour TP'{L) = L -\- 0{L^)^ following from 
the fact that N[E) is a staircase function. 

In figure 6 we show the number variance in the odd symmetry case for L < 50 and different 
values of i^: E = 1000,4000,7000,10000 {6 = 1800). In table I the saturation values S^(i^), 
calculated as the average of E) over L G [15, 100], are listed. It is seen that the agreement 

with the GOE expectation (dashed-dotted line) is restricted to L < 2 ... 7, depending on E. 
Analogous results hold in the even symmetry case. 

It is now well established that the medium- and long-range statistics measured by the 
number variance strongly depend on the (non-universal) short periodic orbits [13, 14, 15, 35, 
16, 17]. Therefore the universality regime, which is expected to be governed by the random- 
matrix model, is restricted to very small correlation lengths Z, which is confirmed by our 
numerical results. In order to enlarge the range of agreement with the random-matrix model, 
the center of the averaging window has to be in the deep semiclassical regime corresponding to 
extremely high lying energy levels. 

According to Berry's semiclassical analysis [14, 15] E^(£') should increase with increasing 
energy E in the limit of large E. A look at figure 6 and table I shows that such an overall 
increase is indeed observed within the relatively large fiuctuations of the number variance. 



9 




Q I I I I I I I I I I I 

5 10 15 20 25 30 35 40 45 50 

L 

Figure 6: Number variance TP'{L,E) for odd symmetry at different energies E: E = 1000 (full 
line), 4000 (long dashed line), 7000 (short dashed Hue), 10000 (dots); 8 = 1800. Dashed-dotted 
line: GOE expectation Yjqq^[L). 



E 


^1{E) (odd) 


E^(i) (even) 


Aoo(i) (odd) 


Aoo(-E') (even) 


1000 


0.59 


0.60 


0.29 


0.29 


4000 


0.75 


0.75 


0.38 


0.37 


7000 


0.79 


0.83 


0.40 


0.41 


10000 


0.85 


0.83 


0.43 


0.38 



Table 1 : Saturation values (£") of the number variance, and Aqo (£") of the rigidity at different 
energies E. Notice that the last digit is uncertain. 

3.5 Spectral Rigidity 

Another statistic measuring two-point correlations is the spectral rigidity [10] 

A3(iv,i^) := (min- / de [N{E + e) - a - be]'^) , (9) 

\ (a.b) L J 



^ -L/2 



E,S 



which is the average mean square deviation of the spectral staircase function from the best 
fitting straight line over an energy range of L mean level spacings. The brackets (. . .)^ ^ again 
represent a local average. 

For short correlation lengths Z ^ 1 the rigidity is independent of the underlying spectrum, 
and obeys ^^{L) ~ j^L^ which is a consequence of N{E) being a staircase function. 

In the theory of random matrices (RMT) TP'{L) and ^^{E) are related through 



J {L'-2L\ + r')j:ly,^ir)dr . (10) 
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By this equation it is possible to obtain Af^^{L) or Af^^^{L) by numerical integration of 
equations (7) or (8). 

In figure 7 the rigidity calculated from equation (9) is shown for different energies E and 
6 = 1800. One clearly observes a saturation effect of the rigidity for large L as it is suggested 
by Berry's analysis [13] of the spectral rigidity using the semiclassical trace formula. This is in 
contrast to the logarithmic increase predicted by random-matrix theory, see equation (7). The 
rigidity is smoother than the number variance, but reaches the saturation later in comparison 

with j:\l,e). 

0.4 I 1 1 -7^ 1 1 1 




' ' ' ' ' 1 

20 40 60 80 100 

L 



Figure 7: Rigidity A.^{L,E) for even symmetry at different energies E: E = 1000 (full line), 
4000 (long dashed), 7000 (short dashed), 10000 (dots), 8 = 1800. The dashed-dotted line shows 
the expectation of random-matrix theory in the GOE case Af^^{E). 

The prediction of the random-matrix model for the rigidity suffers from the same drawback 
as in case of the number variance, i.e., that it is valid for small correlation lengths L only. 

For the following section it is of great importance to know the dependence of the saturation 
value Aco{E) of the rigidity on the energy. Aco{E) is defined as Aco{E) = lim/,^oo A^{L,E). 
According to Berry's semiclassical analysis [13] one expects in the semiclassical limit — )• oo 
for systems with time-reversal symmetry 

Aoo(^) = ^log^ + C . (11) 

The constant C has been estimated by Berry as C = C(/o) = log(47reJ(£')//o) — | where 

d{E) = ^ is the mean level density, /q denotes the length of the shortest periodic orbit and 
e is the base of the natural logarithm. In case of the cardioid billiard we have /q = 2.598 . . . and 
thus obtain C(/o) ~ —0.034 for both symmetry classes. Due to the finite number of available 
energy levels we cannot perform the limit Z — )• oo. Therefore we have determined Aoo{E) from 
a fit of A^{L, E) in the range L G [15, 400] for fixed E to the function 

Af (X) = (^1 + £ + , (12) 
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where Aqo, a and b are fit parameters. We have chosen 9 different energies E from the intervai 
[1000,10000] and determined the corresponding saturation vaiues Aoo(-E') using the fit (12). 
Fitting in turn the asymptotic behaviour (11) to these vaiues we obtained for the constant 
C = —0.048 in the odd, and C = —0.052 in the even case. The resuit of the saturation values 
together with the asymptotic curve (11) are represented for the odd case in figure 8. Notice 
that the values obtained for the constant C deviate from the simple estimate using only the 
shortest periodic orbit. 

0.5 

0.4 

0.3 

0.2 

0.1 



2000 4000 6000 8000 10000 

E 

Figure 8: Saturation values Aoo(-E') (stars) for odd symmetry at different energies, and the 
asymptotic curve (11) with C = —0.048 (full line). 

Another interesting aspect which can be deduced from equation (10) is that there exists a 
simple relation between the two saturation values [35] 

= 2Ao. • (13) 

This relation is confirmed by our numerical data within the expected errors, see table 1, except 
at = 10000 in the even case. 

3.6 The Mode-Fluctuation Distribution P{W) 

As we have seen in the preceding sections, the agreement of the computed number variance and 
rigidity with the predictions of random-matrix theory holds only at small correlation lengths. A 
novel quantity which can be used as an indicator of quantum chaos in spectra is the distribution 
P{W) of the mode fiuctuation W{E) [5, 6]. The function W{E) is defined as the normalized 
fiuctuations of the mode number N[E) around the mean mode number N[E) 

WiE) := ^[i^)^m = . (14) 

^/aZ{e) ^/A^{E) 
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odd symmetry 


even symmetry 


conjecture (15) 




— n nm 




n 


Variance 


0.992 


0.996 


1 


Skewness 


-0.021 


0.008 





Kurtosis 


0.112 


-0.089 






Table 2: First moments of the mode-fluctuation distribution P[W) for the odd and even case. 

By definition Nosc{E) fiuctuates around zero, see figure 3, and furthermore it can be shown [19] 
that the second moment of Nosc{E) tends asymptotically to the saturation value Aoo(-E') of the 
rigidity. 

In [5, 6] the conjecture has been put forward that classically chaotic systems should display 
a universal Gaussian behavior in the limit E ^ oo 

PGaussiW) = ^e-^'^' , (15) 

whereas classically integrable systems should display non-Gaussian distributions P[W). 




-4-3-2-10 1 2 3 4 

W 

Figure 9: Histogram of the mode-fiuctuation distribution P[W) (odd case). Dashed line: 
Gaussian limit distribution. 

Numerically the distribution function was calculated by randomly choosing 10'' energy values 
E G [0, 11000] and determining the histogram of W{E). As one can see in figure 9 the agreement 
between the histogram and the Gaussian normal distribution is very good. In table 2 the mean 
value, variance, skewness and kurtosis for the odd and even cases are shown. 

In order to quantify the visual impression, we applied the Kolmogorov-Smirnov test to the 
cumulative mode-fiuctuation distribution C{W) (see figure 10). We have chosen equidistantly 
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-4 -2 2 4 

W 

Figure 10: Cumulative mode-fluctuation distribution for the even case in comparison with the 
cumulative version of (15). 



11000 points of the unfolded energy interval [0,11000] and determined the maximum value 
D of the absolute difference between the cumulative distribution and the cumulative normal 
distribution. The signiflcance level V of the distance D is in both symmetry cases above 60%, 
which strongly conflrms the conjecture. 



4 The Trace of the Cosine— Modulated Heat Kernel 

It is well known that the Gutzwiller trace formula for bound quantum systems [12, 1] obtained 
from a semiclassical approximation for the density of states is at best conditionally convergent. 
Based on Gutzwillers semiclassical approximation for the level density a generalized trace for- 
mula has been derived by considering the trace TT[h[H^^^)) of a suitable test function h[p) 
of the square root of the Hamiltonian H [36]. The result obtained in the semiclassical limit 
of Planck's constant /i — )• is a periodic-orbit sum rule, which involves absolutely convergent 
sums and integrals only [h = 2m = 1) 

OO OO 7 

^%o= J dph{p)d{p)+ch{o)+j:j: _^ j-,{h{ki,)} . (16) 



7 k=l 



7 



Here the sum on the left-hand side runs over all quantal eigenvalues (not unfolded !) parametrized 
by the discrete momenta {pn = \/En}- The integral on the right-hand side is the so-called "zero- 
length contribution" and is completely determined by the generalized Weyl's law d{p) = '^^J^ ^ . 
C is the constant term in Weyl's law. 

The 7-summation runs over all primitive periodic orbits 7 of length /-y, whereas counts 
multiple traversals corresponding to periodic orbits of length kl^; is the monodromy matrix. 
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and J-'^{h[x)} denotes the Fourier transform of the test function h[p) (incorporating the phase 
shift due to the Maslov index) 



J^^{h{x)} 



1 



dph[p) cos[px — kp^^ 



:i7) 



where denotes the Maslov index which is twice the number of reflections at Dirichlet bound- 
aries plus the maximal number of conjugate points fi^ (see e.g. [37, 38]). 

Starting from the generalized trace formula (16) we can study the inverse problem of quan- 
tum chaos: to obtain the lengths of the periodic orbits from the quantal energy levels. This 
problem can be solved by considering the trace of the cosine-modulated heat kernel [39, 40] 
Tr|cos (j^—Ay^^L^ which is obtained from equation (16) by choosing the following test 

function 

h{p) = cos{pL) e-P'' ; p = Ve , t>0 . (18) 

Due to the geometry of the cardioid billiard with its focusing boundary, the orbits listed in 
table 3 have a non-vanishing maximal number of conjugate points fi^. Thus can be odd or 
even which results in a sine- or cosine- Fourier transform respectively. For even kp^ = 2m one 
has a cosine- Fourier transform which yields for (18) 



J^^{h{x)} 



L] 



exp 



4t 



-\- exp 



[x + L] 



(19) 



whereas for odd kp^^ = 2m -\- 1 the result of the sine-Fourier transform of (18) is given by 



J^^{h{x)} 



1 

Vt 



L] 



it 



+ (x + Z)iFi 1, 



jx + L) 



(20) 

where \Fi{a, 6, z) is Kummer's function. 

The Fourier transforms (19) and (20) clearly display the different kind of shapes of the 
periodic-orbit contributions to the trace formula as a function of L for fixed t. If kp^^ is even, 
one obtains a Gaussian peak at Z = kL^, whereas for odd kp^^, one has a zero at the lengths of 
the periodic orbits because 1-^1(1, |,0) = 1. This is similar to the observation made by Sieber 
et al. [41] in case of the Fourier analysis of the spectrum. 

The trace of the cosine-modulated heat kernel allows one to extract the length and stability 

of the first orbits and even the value of kp^^ modulo 4. The last property cannot be read off in 

2 

case of the often considered so-called power spectrum D{x) = /q™^'' dp e^P^[d{p) — d{p)] . 

With the first 11000 levels and a smoothing parameter t = 0.0001 we obtained the graph 
(full line) shown in figure 11 for the odd symmetry case. For a proper comparison with the 
trace formula, we have to subtract from the right-hand side of (16) the contribution of the 
unknown part of the energy spectrum, {En, n > N = 11000}. This contribution is estimated 
by first replacing the sum over the discrete spectrum by a Riemann-Stieltjes integral, and then 
approximating the integral by means of Weyl's law, equations (3) or (4): 

00 °;? 7 - 

h{pn)= / dN{E)h{^/E) ^ / dph{p)d{p) . 



n=N+l 



Subtracting this term from the right-hand side of (16) implies that the integration in the 
integral in (16) has to be restricted to the interval [0, \/En]. With this approximation we also 
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Figure 11: Trace of the cosine-modulated heat kernel for odd symmetry with t = 0.0001. Full 
line: sum over the first 11000 energy levels, dashed line: periodic-orbit sum. The squares 
indicate the lengths of the first primitive periodic orbits, see table 3 and figure 12, and of 
their repetitions. 



show in figure 11 the right-hand side of (16) (dashed line) by using the geometrical data of 
the primitive periodic orbits known so far (see table 3, and figure 12). We observe a striking 
agreement between both graphs which indicates the accuracy of the computed energy levels. 
Notice the excellent resolution of the two peaks near L 6.6 which correspond to orbits No 4 
and 5. 

In the interval L G [7.1,7.6] there are some differences between the trace of the cosine- 
modulated heat kernel calculated from the eigenvalues and the periodic orbit sum respectively, 
which are probably caused by orbits which run into the cusp of the cardioid billiard. Their 
contribution is expected to be of higher order in h and thus is not included in the trace formula 
(16). Similar differences exist in the even symmetry case. A detailed analysis of the trace 
formula for the cardioid billiard is in progress and will be presented in a separate publication. 
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Orbit No 


7 


Tr AL, 




//^ (D) 


^1 (N) 


1 


2.598 


-2.50 


1 


5 


3 


2 


4.618 


-4.40 


2 


8 


6 


3 


5.918 


-9.74 


3 


11 


9 


4 


6.585 


7.62 


3 


13 


9 





D.D / 


-ZO.OD 


/I 

4 


1 A 
14 




6 


7.103 


12.45 


4 


16 


12 


7 


9.552 


-16.31 


4 


18 


12 



Table 3: Dynamical data of the first primitive periodic orbits. is the length of the orbit 7, 
Tr Af^ denotes the trace of the monodromy matrix, is the maximal number of conjugate 
points, and z/^ is the Maslov index (for odd (D) and even (N) symmetry). 
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5 Summary and Discussion 



In this paper we have presented a detailed analysis of the spectral statistics in the cardioid bil- 
liard. This dynamical system is strongly chaotic and thus constitutes an ideal testing ground 
for quantum chaology. It is found that the level-spacing distribution is in very good agreement 
with the Wigner distribution. However the number variance and rigidity clearly violate the 
universal laws of random-matrix theory for medium and large correlation lengths. The ob- 
served non-universal behaviour agrees with semiclassical expectations. Our results thus further 
strengthen the evidence accumulated so far showing that the applicability of the random-matrix 
model is limited to short-range correlations only. 

A promising candidate for a universal measure of quantum chaos in spectra is the mode- 
fluctuation distribution P[W). We have shown in sec. 3.6 that this novel spectral statistic 
indeed displays at a high significance level a universal Gaussian behaviour in agreement with a 
recent conjecture. In view of this result it appears even more urgent to inquire for a rigorous 
derivation of the conjecture from first principles. 

Finally, in section 4, we have studied the trace of the cosine-modulated heat kernel and have 
compared it with the theoretical result derived from Gutzwiller's semiclassical trace formula. 
The overall agreement between "experiment" and theory is striking. It is amazing to see in 
figure 11 for instance the pronounced zero at the location of the length of the shortest periodic 
orbit. The observed structure follows precisely the theoretical curve as described by equation 
(20). This is a nice illustration of "inverse quantum chaology" : Knowing the quantal eigenvalues 
we are able to determine not only the lengths of classical periodic orbits, but also whether they 
possess conjugate points or not. This agreement with periodic-orbit theory calls for a more 
elaborate investigation of the trace formula in the case of the quantized cardioid billiard. Work 
along this line is in progress and will be published later. 
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